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Abstract. Using a reasonably large chemical reaction 
network consisting of 421 species, we show that along with 
normal chemical evolution of molecular cloud during col- 
lapse and star formation, significant amount of adenine, a 
DNA base, may be produced after an evolution of 10^^^ 
years. This findings may shed light on whether life on earth 
had to begin from scratch or these molecules could have 
contaminated the earth from the beginning. 
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1. Introduction 

Understanding origin of life on earth is a very challenging 
astrophysical problem of today. Study of isotopic ratios 
on carbon inclusions in ferric salt deposits in Greenland 
shows that life existed on earth before 3.85 billion years 
(Mojzsis at al. 1996; Holland 1997; Schidlowski 1993) or 
possibly even before 4 billion years (Mukhin & Gerasimov 
1993). On the other hand, data from isotope systematics 
show that the earth is between 4.55 to 4.57 billion years 
old (Faure 1986). The time available for life to develop 
on the earth from scratch is possibly around 0.5 — 0.6 
billion years which is very short (Crick & Orgel 1973). 
One of the ways the formation of life on earth could be 
assisted if the building blocks of life, namely, amino acids 
and at least some of the bases of DNA could be produced 
well before the planetary formation, presumably during 
chemical evolution of molecular clouds. Whereas extensive 
work is present in the literature to study evolution of some 
of the organic and inorganic molecules, no work is reported 
so far which study the formation of bio-molecules or pre- 
biotic molecules such as amino acids, sugars, DNA, RNA 
etc. 

Bases of DNA such as adenine can be pro- 

duced by successive addition of HCN molecules in four 



steps (Volkenshtein, 1983). In the present letter we per- 
form a hydrodynamic simulation and followed chemical 
evolution during the molecular cloud collapse and show 
that a significant amount of adenine (with mass fraction 
^adenine ~ 10~^°) could bc produccd. This indicates that 
planets and comets could have formed with contamina- 
tions of pre-biotic molecules. If correct, this may resolve a 
long-standing problem of origin of lives on earth. 

2. Reaction Netvifork 

We take the UMIST database (Millar, Farquhar & Willacy 
1997; hereafter referred to as MFW97) as our basis of 
chemical reactants and reactions, but added several new 
reactions such as synthesis of amino acids (alanine and 
glycine), hydroxy-acids (glycolic and lactic acids), DNA 
base (adenine, see, Volkenshtein 1983), urea synthesis etc. 
These new reactions make the total number of species to 
be 421. The rate co-efficients of these additional reactions 
are difficult to find, especially in the environs of a molec- 
ular cloud. To use UMIST database, the rate constant for 
a two body reaction is written as (MFW97), 



k(T/300)^exp(-7/r) 



cm'^ s 



(1) 



where, a, (3 and 7 are constants and T is the tempera- 
ture. Amino acid synthesis rate was estimated from Fig. 
8 of Schulte & Shock (1995). Urea synthesis rate is kept 
comparable to the rates given in UMIST table. The rate 
constants were taken to be a = lO""'^", /9 = 7 = for 
each two-body reactions. The rate constants for adenine 
synthesis was chosen to be similar to other two body re- 
actions [a — 10^^" /? = 7 = for each HCN addition in 
the chain HCN ^ CH{NH)CN NH2CH{CN)2 
NH2{CN)C = C{CN)NH2 ~* H5C5N5 (adenine)]. This 
should provide a sufficiently good estimate for the final 
abundance of these pre-biotic molecules. 

Initial composition of the cloud before the simulation 
begins is kept to be the same as in MFW97, and for- 
mation of H2 is included using the grain-surface reac- 
tion with rates as in MFW97. The initial mass fractions 
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are taken to be the same as in MFW97 (but converted 
to mass fractions), i.e., H:He:C:N:0:Na:Mg:Si:P:S:Cl:Fc 
= 0.64:0.35897:5.6 x 10-4:1.9 x 10-'*:1.81 x 10-^:2.96 x 
10"^: 4.63 X 10-^:5.4 x 10-^:5.79 x 10-^:4.12 x 10-^:9 x — 
10-^:1.08 x 10-8. H 



3. Hydrodynamic considerations 

We use simple molecular cloud model witti generic ini- 
tial conditions: initial size of the molecular cloud ro = 
3 x lO^^cm, average temperature of the cloud T = lOK, 
and angular velocity of the cloud = 10"^^ rad s"^. The 
sound speed becomes = (fcT//imjj)^/^ ~ 19200cm 
and corresponding initial density (Shu, Adams & Lizano 
1987) is /9 = a^/27rGr^ = 10-^^g cm"^ and accretion rate 
is M = 1.06 X 10^°g s~^. Here, jjL = 2.3 mean molecular 
weight, m,H is the mass of the hydrogen atom, and k is the 
Boltzmann constant. In the isothermal phase of the cloud 
collapse, density p oc (Chandrasekhar 1939) and the 
velocity is constant. When opacity becomes high enough 
to trap radiations (say, at r = rtr), the cloud collapses 
adiabatically with p oc r~^/^. In presence of rotation, cen- 
trifugal barrier forms at r = Tc, where centrifugal force 
balances gravity. Density falls off as p oc r~^/^ in this re- 
gion (Hartmann, 1998). Following Shu, Adams & Lizano 
(1987), we compute the density, temperature and velocity 
distribution inside the cloud and follow the chemical evo- 
lution at the same time. With density chosen as above, 
the initial constant velocity of infall becomes 8900cm 
and below r = rg velocity oc r^^/-^ to preserve the accre- 
tion rate in a disk like structure of constant height. Since 
for the parameters chosen (generic as they are) Tc > rtr, 
we chose T oc 1/r inside the centrifugal barrier (r < Tc) 
as in an adiabatic flow. We follow the collapse till a ra- 
dius of lO^^cm is reached. Figure 1 shows the temperature 
and mass density (in units of lO^'') distribution in a log- 
log scale. Our choice of O is small by at least an order 
of magnitude from that of the observed value in some of 
the clouds. This was done to obtain atmosphere at a dis- 
tance of around lAU from the proto-star to be as close to 
reducing {C/O > 1). 
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Fig. 1 : Log-log plot of temperature (dotted) and mass density 
(in units of 10'^'') as functions of the radial distance. Vertical 
dashed line is drawn at lAU. 
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Fig. 2a : Log-log plot of the mass fractions of some of the 
lighter and intermediate mass species as functions of the ra- 
dial distance. Alternate species have been plotted with dotted 
curves for clarity. Upper axis shows time elapsed in seconds 
since collapse began. Vertical dashed line is drawn at lAU. 



4. Results 

Figure 2a shows the log-log variation of mass fractions of 
some of the light and intermediate species of our simula- 
tion with radial distance (in cm). For clarity, we plot the 
curves alternately by solid and dotted type. In Fig. 2b, 
we plot the variation of more complex molecules. Since 
they start with zero abundances, we plot them from ra- 
dius 2.5 X lO^^cm for clarity. We put a long dashed vertical 
line at lAU, the distance of the earth with respect to the 
sun. On the upper axis, we have put time (in seconds) 
elapsed since the beginning of collapse at radial distances 
of lO^^cm, lO^^cm, lO^^cm and lO^^cm respectively. To- 
wards the end of the collapse, time spent is negligible. 



Because of low initial velocity, H, N, C and O are de- 
pleted much before lOOAU and HCN, NH3, CH3CN etc. 
rises as also the more complex molecules (Fig. 2b) which 
form out of them. Inside Tc, as matter falls faster and 
spends lesser time, the depletion of lighter molecules are 
controlled, until the density and temperature also rises so 
high that depletion started once more. Inside. 3 x lO'^^ cm, 
the composition changes in a very short timescale Around 
r ^ lOAU, the mass fractions of adenine, urea and glycine 
arc already significant. At r = lAU, Xadenine = 6 x 10"'^"'^. 
Since the mass of the earth is around 4 x lO^^g, this cor- 
responds to 2.4 X lO^^g of adenine which could have con- 
taminated the earth at the time of formation (it could be 




Fig. 2b : Log-log plot of the mass fractions of a few complex 
molecules (marked) as functions of the radial distance. Upper 
axis shows time elapsed in seconds since collapse began. Verti- 
cal dashed line is drawn at lAU. 

higher since metallie content of earth is mnch above the 
average molecular value. Also, one has to include the effect 
of dust-gas mixture in the molecular cloud) This compu- 
tation, does not consider the destructions of adenine at 
a higher temperature region, and it is likely that much 
of these contaminants are destroyed during collapse and 
formation of proto-carth. However, comets formed in the 
inner cloud could carry away these pre-biotic molecules 
and deposit them during future impacts on planets. It is 
to be noted that around lAU, the composition is close to 
reducing (C/O > 1) type (Fig. 2a) which is favourable for 
the formation of bio- molecules. 

Since regions of low density of molecular cloud could 
have a very long evolution time scale, it may of interest 
whether pre-biotic molecules could have formed in a very 
low density isothermal region. To this effect, we choose a 
cloud of mass density p = 1.6 x lO^^^g cm~^ and T = lOK 
and let it evolve for 10^ years. Figure 3 shows the results of 
this simulation where we plot the variation of abundance 
with time (in seconds). We find that the abundance of 
adenine, for instance, is around 5 x 10~^^. This is a very 
generic condition, and the abundance is significant. We 
therefore believe that adenine could be produced during 
the molecular cloud collapse. 

5. Concluding remarks 

It is generally thought that molecules vital to living sys- 
tems on earth may have started from scratch. In this Let- 
ter, we have explored the possibility of formation of a DNA 
base, namely adenine, and a few of the simpler amino acids 
during a generic molecular cloud collapse and show that 
mass fraction of adenine is significant which might have 
contaminated the planetary disk even before the forma- 
tion of the planets themselves. In a static evolution of a 
cloud of constant density and temperature, the formation 
of adenine is also important, and therefore our conclusion 
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Fig. 3: Log-log plot of the mass fractions of some pre-biotic 

molecules (marked) as functions of time during the evolution 
of a static cloud of constant density and temperature for 10® 
years. 

may be; very generic. The rates wc; have chosen arc reason- 
able, though uncertainties of an order of magnitude cannot 
be ruled out at this stage as more accurate estimates re- 
quire one to have temperature dependent rate coeflacients 
of both forward and reverse reactions. 

Since earth is a generic planet and since it seems ap- 
parent that bio-molecules are as naturally produced as 
any other complex molecules, we believe that DNA bases 
produced in collapsing cloud could have contaminated the 
earth, and at the same time there should be many such 
planets in each galaxy where DNA based lives should 
flourish. It is not yet clear whether or not these molecules 
survived during further evolution of planetary disks lead- 
ing to eventual formation of planets or actual contami- 
nations were made much later by comets carrying away 
protecting these molecules. These vital issues will be ad- 
dressed in due course. 
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